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(54) Color seperation method and apparatus for same 

(57) A colorant selection mechanism for a color 
reproduction device is disclosed that is capable to find 
all the sets of colorants that yield a specific color and to 
select the most appropriate one. Both the method and 
an apparatus are disclosed. 
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Description 

Field of the invention 

5 The present invention relates to devices and methods for image rendering by means of color reproduction devices. 

Background of the invention. 

The independent values with which the color device can be addressed are called colorants or inks. For the purpose 
10 of ger^ality. the values for these colorants can always be scaled to range from 0 to 100 %. A color reproduction device 
with n colorants is called a printer or an n-tnk process. 

With cokxant space rs meant an n-dimensional space with n the number of independent variables with which the printer 
can t)e addressed In the case of an offset printing press, for example, the dimension of the, space corresponds with 
the nurrtoei of mKs of the printer 

15 With cdof «pace is meant a space that represents a number of quantities of an object that characterize its cobr. In most 
practical situatKyis. colors will be represented in a 3-dimensional space such as the CIE XYZ space. However, also 
other characteristics can be used such as multispectral values based on filters that are not necessarily based on a lin- 
ear transJormabon oi the color matching functions to represent color. A typical example is an n-dimensionat space of 
which the ames correspond to densities. 

20 With colorant garrwl or colorant domain is meant the delimited region in colorant space of colorant combinations that 
are physicaffy rM«zabi« by a given printer, taking into account possible extra limitations on colorant combinations. In 
this invention i « s^^jposed ttiat colorant limitations can preferentially be specified by mathematical formulas that 
describe the boundaries ot the desired colorant gamut. 

The colorant garr^ of a crryk offset printing press for example is often delimited by a linear condition tiiat limits the sum 

25 of the four cotar arte (tor •sample to 340%). 

A printer model « m mathematical relation that expresses color values as a function of colorants for a given printer. The 

variables of tt>e cctorams are denoted as Ci . C2 Cn with n being the dimension of the colorant space. 

A printer or nnnk process « completely characterized by its colorant gamut witii a number of colorant limitations and tiie 
printer model By rvertvig the printer or the n-ink process is meant that the printer model is inverted. If an n-ink process 

30 is given with cokyani imtations in the colorant domain and if an m-ink process is deduced from this n-ink process by 
setting n-m colorams 10 constant values In the printer model of the n-ink process, the colorant limitations are inherited 
by this nvink process d the cokKant gamut of the m-ink process is restricted by colorant limitations obtained by setting 
the n-m colorants w>e colorant limitations of the colorant gamut of the n-ink process to their constant values. 
Extracting an nvr* process out of an n-ink process with m < n. means that n-m out of the n colorants in the printer 

35 model are replaced a constant value. The colorants of such an m-ink process vary between the minimum and max- 
imum values they can have according to the colorant gamut of the n-ink process and the colorant limitations are not 
inherited t)y the aytnk process The m-ink process is called the extracted m-ink process of the n-ink process. 
If an m-ink process is extracted from an n-ink process, and the n-m out of the n colorants in the printer model of the n- 
ink process are replaced by their minimum or maximum value as defined by the colorant gamut of tiie n-ink process, 

40 the nvink process is caMed an extracted boundary m-ink process of the n-ink process. In general there are 

{n-m) \m\ 

45 

extracted kxHirxtery rTv«n*i processes of an n-ink process. The colorants in the colorant gamut of an extracted boundary 
m-ink process vary betwroen the minimum and maximum values they can have according to tiie colorant gamut of the 
n-ink process and the colorant limitations are not inherited by the m-ink process. With color gamut is meant a delimited 
region in color space, containing colors that are physically realizable by a given printer, taking into account possible 
so extra colorant limitations Typical devices we have in mind are offset or gravure printing presses, thermal dye sublima- 
tion processes, therrr^ wax transfer processes, electrophotographic processes and many otiier processes. Neverthe- 
less, tiie invention can also t>e applied to other color reproduction devices such as color displays, color photography or 
film recorder 

A typical printer is capable of laying down various amounts of cyan, magenta and yellow colorants. The modulation of 
55 these various arTX)ur<s can be done by altering the concentration of the dye that is laid down onto tiie substrate or by 
means of one of the haiftorwig technques. an overview of which is found in the article by J.C. Stoffel and J.F. Moreland 
"A survey of Electrons Tec^m^iues for Pictorial Reproduction". IEEE Transactions on Communications. Vol. COM-29. 
No. 12. Dec. I98t More and updated information on halftoning is found in the patent literature and a list of recom- 
mended patents is tound m the artide by Peter R. Jones, "Evolution of halftoning technology in the United States patent 
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literature". Journal of Electronic Imaging 3(3), 257-275. July 1994. 

The calculation of the correct amounts of colorant for the rendering of color images on a printer is called the color sep- 
aration problem. Most of the color separation strategies known In the art comprise the following steps: 
In a first step, the relation between the amounts of colorants and the resulting colors, on a printer is characterized. This 
£ is done by first printing a set of colorant combinations that spans the dynamic range of the printer and measuring the 
resulting colors. An example of such a set is the ANSI IT8.7/3 form. 

In a second step this relation is mathematically modeled, to obtain the printer model. The printer model consists of 
some form of an analytical expression that predicts color for a given combination of colorant amounts. A well known 
mathematical model is the set of "Neugebauer equations". 
w In a third step the printer model is Inverted. This is necessary as the color separation problem Is involved with finding a 
set of colorants that renders a given color and not vice versa. Model inversion is usually done by means of some itera- 
Ifve search strategy, such as the Newton Raphson iteration. 

A difficulty is that many colors that can occur in photographic images are not reproducible at all by a given printer and 
hence that no physically realizable set of colorants can be found to render these colors. The solution to this problem 
15 consists of first converting the "non reproducible color" Into a "reproducible color". Up till now there exists no universally 
accepted strategy to achieve this goal, and the evaluation of any such strategy is therefore necessarily based on aes- 
thetic and therefore subjective judgment A discussion of strategies to deal with the problem is found in the article by 
Jim Gordon "On the rendition of unprintable colors", 1987 Proceedings of the Technical Association of the Graphic Arts 
pp 186-195 

:v Because it will help to understand the problems in the prior art. a few more words of explanation are given on the Neu- 
gebauer equations for a 3-ink process and a few methods are presented that are currently in use for their inversion. The 
generatizaton of the Neugebauer process for n inks is straightforward. 
The Neugebauer equations for a 3-ink process. 

The printing with three inks and three halftone screens results theoretically in 8 possible combinations of ink overlap, 
2s The Neugebauer expressions predict the resulting color as a linear function of the colors of these combinations. The 
N ugebauer equation for the X tristimulus value in a 3-ink process is: 

X(C,.C2.C3)=a^X^+aiX,+a2X2+a3X3+a23X23+ai3X;3+ai2>^i2+ai23^i23 

30 The terms X,^ in the above equation are the X tristimulus values of the corresponding overprints. The Neugebauer 
expression tor the Y and 2 tristimulus values are obtained by replacing the X tristimulus values by the corresponding Y 
and Z values respectively. If it is assumed that the relative positions of the halftone dots is random, the Neugebauer 
coeff icients can be calculated from the Demichel equations that predict the fraction of each combination of the three 
inks as a function of their respective dot percentages Ci, C2 and C3 : 
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aw=0-Ci)(l-C2)(1-C3) 
ai=(Ci)(1-C2)(l-C3) 
a2=(i-Ci)(C2)(i-C3) 
a3=0-Ci)(l-C2)(c3) 

a23=(1-Ci)(C2)(C3) 
ai3=(Ci)(1-C2)(C3) 
ai2=(Ci)(C2)(1-C3) 

Substitution of the Demichel equations in the Neugebauer equations and rearranging the terms leads to the following 
representation of the three Neugebauer equations for a 3-ink process: 

X = ko+kiCi+k2C2+k3C3+ki2CiC2+ki3CiC3-fk23C2C3+ki23CiC2C3 
Y = lo+liC,+l2C2+l3C3+l 12C 1 C 2+1 13C 1 C 3+1 23C2C 3+1 123^,0203 
Ys= mo+m,Ci+m2C2+m3C3+mi2CiC2+mi3CiC3+m23<^2<53+"^123CiC2C3 
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These equations predict the XYZ tristimulus values as a function of the amounts of the three colorants c^, and C3. 

Interpretation of the Neugebauer equatbns 

5 It is worthwhile to consider a number of different mathematical interpretations of the Neugebauer equations. 

According to a first interpretation, the Neugebauer equations can be seen as a set of three (one for X, Y and 2) trilinear 
interpolation formulas. This becomes immediately clear by looking at the form of the coefficients of the Neugebauer 
equations. These coefficients result in exact predictions of X, Y and Z for values of c^. Cg and C3 when they are 0 or 
100%, and in trilinearly interpolated values for other values. If the Neugebauer equations are used as multilinear inter- 
to polation formulas, they can be used to model any color de\nce in any color space. 

Instead of using coefficients that make the Neugebauer equations exact for 8 colors, it is sometimes desirable to have 
coefficients that minimize the error (for example in terms of "least mean square") over a large set of colors. In that case 
a regression technique can be used to obtain the Neugebauer coefficients. The equations in that case have the nature 
of an "interpolating polynomial" . 
15 An interesting observation is that for a given X. Y and Z, each of the three Neugebauer equations represents a "surface" 
in the 3-dimensional c^ , C2, C3 space. As Figure 1 shows, can the solution of the Neugebauer equations for a given X, 
Y and Z value be seen as the intersection between the three surfaces that correspond to these values. 

Improving the precision of the Neugebauer equations 

20 

Real printing processes seldom behave exactly according to the physical model on which the Neugebauer equations 
are based, and this explains the deviations that can occur between the predicted and measured color for a given color- 
ant combination. 

The interpretation of the Neugebauer equations as an interpolating polynomial leads to the introduction of additional 
25 higher order terms to improve their precision. In theory it is possible to increase the number of terms until any desired 
precision (in terms of "least mean square error") is obtained, but practical arguments such as numerical stability and 
the amount of computation needed for the calculation of large sets of coefficients from large sets of data put constraints 
on the achievable precision. 

An alternative approach, and this is our prefen-ed one, is to "localize" the coefficients of the Neugebauer printer model, 
30 and this leads to a "piece-wise" interpolation approach known under the name of the localized or cellular Neugebauer 
equations (LNE). 

The different subdomains are called the Neugebauer cells. In Figure 2 the colorant domain of a cmy process is pre- 
sented that is divided into 8 Neugebauer cells. Per Neugebauer cell, the printer is modeled with a 3-ink Neugebauer 
process. 

35 Various other modifications of the Neugebauer model have also been suggested to improve their precision, among 
which the n-modif led and the spectral Neugebauer equations are dted. These and other models are discussed in the 
overview article by R. Rolleston and R. Balasubramanian "Accuracy of various Types of Neugebauer Model", presented 
at the 1993 IS&T and SID's Color Imaging Conference on Transformation & Transportability of Color and in the article 
by Kang H. R.. "Applications of color mixing models to electronic printing". Journal of Electronic Imaging 3(3), July 1 994. 

40 

Inversion of Neugebauer equations by means of iterative methods 

The Neugebauer equations predict color from colorant. Separation is involved with the opposite transformation. Several 
methods have been known to achieve this "model inversion", and all of them are based on an iterative search that con- 
45 verges to the desired solution. One example is the Newton Raphson inversion method. Since the derivative of the Neu- 
gebauer equations exists and is easily calculated. Newton Rapfison iteration can be used. The values off Cf. C2 and C3 
in the (i+1)'^ iteration are found from the values in the i^ iteration as follows: 
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55 in which the corrections Aci , AC2 and AC3 are found by solving the set of linear equations: 
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dX/dc^ dX/dc^ dX/dc^ 
dY/dc^ dY/dc^ dY/dc^ 
dZ/dc^ dZ/dc^ dZ/dc^ 

10 

The values. AX. aY and AZ represent the differences between the XYZ values of the desired color and the color that is 
predicted by the Neugetauer model for c^. C2 and C3 in the ith iteration. The iteration is stopped' when the values of . 
C2 and C3 have become stable. 

The Newtonian iteration converges quickly to the desired solution provided that a reasonable start value was selected. 

15 Simple as this sounds, however, this is not always obvious, and several complications may occur. 

A first problem relates to the possibility that at some point in the iteration at least one of the partial derivatives of the 
color components becomes zero or very small. A typical example of this occurs in the cmy process in that part of the 
colorant space where each of the three colorant components is near its maximum value (100%): the Y color component 
varies very little for changes Ac. Am and Ay on the amounts of colorants c. m or y in that case. A situation like this can 

20 bring the set of equations to near singularity, and can result in overly large corrections Ac. Am and Ay for the colorants, 
causing the iteration to diverge. This is not unsimilar to the example of the Newton Raphson iteration on a 1 -dimen- 
sional function in Figure 3. 

Lesser known and certainly lesser understood are the problems that can occur with the model inversion of a 3-ink proc- 
ess where the same color can be rendered with more than one set of colorants. As Figure 4 shows, the graphical inter- 
ns pretation of this situation is that the three surfaces that correspond to each of the Neugebauer equations have more 
than one intersecting point in the colorant gamut. The drawing was based on the behavior of a printing process with 
cyan, yellow and green inks on a white substrate. It is true that the existence of multiple solutions, that fail within the 
colorant gamut, is relatively uncommon for "well-behaved" printing processes such as the cmy offset printing process, 
and this explains probably why the problem has not got much attention up till now. Our experience shows, however, that 
30 multiple solutions are far less exceptional when unconventional processes are to be modeled. Examples include a ther- 
mal dye sublimation process where above certain dye concentrations in the substrate an "inverse" flow of dye occurs, 
leading to non-monoticity and double solutions for certain colors. Another example is the printing process used in many 
commercial telephone directories where cyan, magenta and black inks ai-e used to print on a yellow substrate. It can be 
readily understood, even intuitively, that more than one combination of the three colorants can be found that yields cer- 
35 tain colors. Yet another example is found in the field of the "HiFi color printing" processes where, besides cyan, 
magenta, yellow and black inks, in addition red, green and blue inks are used: the 3-ink process where, for example! 
cyan, yellow and green inks are printed in the presence of magenta ink leads to "multiple solutions" for certain colors. 
Model inversion of the atwve processes by means of iterations woukJ lead to different sets of colorants depending on 
which "start vector" is used. This effect becomes problematic as images with continuous color transitions can be sepa- 
40 rated into colorant transitions that are discontinuous. At the transitions where the separations "jump" from one solution 
to the other, even small register variations would manifest these discontinuities as disturbing artifacts. 
Opposite to the case of having multiple solutions is the case that the inversion of a set of Neugebauer equations has 
no (real) solution at all. Indeed, as intuition tells and as a more rigorous mathematical analysis would confirm, not alt of 
the XYZ color space can be "addressed" from the c^. Cg and C3 colorant space, even if the colorant values are allowed 
45 to vary from minus to plus infinity. Inversely is it possible that an XYZ color is offered to the iteration process for which 
no (real) set of colorants exists, and hence will an iteration process for such an XYZ cotor never converge. The graphical 
interpretation of this case is that the three surfaces that correspond to each of the Neugebauer equations have no inter- 
section at all. ft may look senseless to try to separate an XYZ color for which no set of colorants exists, but the truth is 
that this is not always known in advance. Many methods, for example, that are used to calculate the reproducible color 
50 gamut of a printing process try to separate a large set of XYZ colors, including colors Ibr which no set of colorants exists. 
All the above problems become only worse when higher order models than the "simple Neugebauer equations" are 
used for the sake of better precision: the additional inflection points that are introduced by the higher order temis only 
increase the chances of multiple solutions or divergence caused by local extremes. 

Depending on the start vector one or the other solution is found or, in the worst case, no solution at all. Never is any 
55 information available on how many other solutions actually exist what they are and which one should be retained. 
Also with other inversion techniques such as simplex methods the same problems occur as with the inversion by means 
of Newtonian iteration, since they are related more to the printer behavior than to the inversion method. Local minimums 
can cause the minimization strategy to converge to the "wrong" set of colorants, and multiple solutions are usually not 
detected. No meaningful results are obtained if the model inversion is done for an XYZ color that has no colorant solu- 
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lion. 

Calculation of separation look up tables. 

The separation of a color into a set of colorants by means of any of the above techniques is a computationally inten- 
sive task When large numbers of colors are to be separated, as is the case for the separation of high resolution images, 
look up tables and interpolation techniques are used. The look up tables have a 3-dimensional address space of which 
the axes are related to the color coordinates of the image. With every addressable position in the table corresponds a 
set of colorants that are calculated off line by means off one of the above strategies. More information of the use of look 
up tables for this purpose is found in the patents US 3,893,166. US 4.334.240. US 4.751 .535. EP 0.550.212 and DE 42 
34 985. and in the article by Kanamori "Color Correction Technique for Hard Copies by 4-Neighbors Interpolation 
Method". Journal of Imaging Science and Technology. Vol. 36, No. 1 , Jan./Feb.. 
A special problem occurs in these tables when two neighboring addressatrfe positions are filled in with colorant values 
that correspond to discontinuous solutions of the Neugebauer equations. 
15 A ver y lar ge en^or on the printed color can originate when interpolation is used to obtain the colorant values for a color 
tfiat lies m between these addressable positions. 

Otjjects of the invention. 

^ H should be dear from the above explanation that the inversion of the Neugebauer equations by means of the existing 
techn^jues does not always converge to a representative set of colorants, and that this can lead to discontinuities and 
color error when images are separated or look up tables are calculated for that purpose. 

The problems are particularly manifest in the case of separations for "unconventional printing processes". Therefore a 
need exists for a more robust solution that, under all conditions, finds all the colorant sets that render a given color and 
25 then retains one of them based on some criterion. Such a method is the object of our invention. More specifically. 

• It IS an obiect of the invention to find the number of possible sets of colorants for the rendering of a color by means 
of an n-ink process in which n is equal to or lower than the dimension of the color space. 

It IS another object of the invention to find these combinations of colorants. 
30 • It IS another object of the invention to select at least one of the combinations of colorants that yields the specified 
cotor. 

• It is another object of the invention to select a set of colorant combinations by minimizing a cost function if the 
dimension of the colorant space is larger than the dimension off the color space. 

• It is another object of the invention to provide an improved separation technology that avoids discontinuities in sep- 
35 arated color images 

• It is another object off the invention to find the colorant set for the rendering of a given color that yields the smallest 
change in color in the presence of variations on the amounts of colorant. 

• It IS another object off the invention to provide a robust method for the inversion of the Neugebauer model for up to 
three colorants in a 3-dimensional color space. 

^ • It IS another object of the invention to provide a robust method for the inversion of the localized Neugebauer equa- 
tions for up to three colorants in a 3-dimensional color space. 

• It is another object of the invention to provide a robust method for the inversion of the Neugebauer model for n col- 
orants (n>s:3) in a 3-dimensional color space. 

• It is another object of the invention to provide a robust method for the inversion of the localized Neugebauer equa- 
45 tions for n colorants {n>=3) in a 3-dimensional color space. 

Further objects and advantages off the invention will become apparent from the description hereinafter. 
Summary off the invention 



so 



The above mentioned objects are realised by the specific features according to claim 1 . Preferred embodiments of 
the invention are disclosed in the dependent claims. 

In the broadest sense. Neugebauer cells are solids, as a portion from a colorant space, such that each point belonging 
to the colorant gamut belongs to at least one Neugebauer cell. Such Neugebauer cells may thus have common points 
55 regions or solids. 

In a specific embodiment, the Neugebauer cells partition the colorant gamut in mutually disjunctive portions. This 
means that the intersection between each two Neugebauer cells is empty or comprises at most one surface. In a more 
preferred embodiment, each Neugebauer cell is a cuboid, i.e. a rectangular parallelepiped. A specific embodiment off a 
cuboid IS a cube. The colorant gamut may be split cubes, all having the same size. 
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In the broadest sense, "Neugebauer equations" as referred to in the current invention are equations having at least one 
non-linear term. A non-linear term is e.g. a quadratic term, cubic, etc. or a cross-product between two or more linear 
variables (e.g. c^ All sets of equations, all having constant or linear terms only, are excluded from the concept of 

"Neugebauer equations". 

5 Preferentially, Neugebauer equations are polynomials that express color values in terms of bilinear, trilinear or quadratic 
colorant values. In a specific embodiment, Neugebauer equations are limited to the formulae for X, Y and Z given under 
the title "The Neugebauer equations" in the background of the invention, which are referred to as trilinear equations. 
One of the key elements 0I our invention is a new and robust method for finding all the colorant sets, including both the 
"rear and the "complex" solutions, that render a given color specified in a 3-dimensional space for a process that is 

70 modeled with the Neugebauer equations for three colorants. The "corrplex" solutions have obviously no physical Inter- 
pretation, neither do these real solutions In which at least one colorant falls outside the colorant gamut. These solutions 
are therefore omitted and need no further consideration. The number of remaining solutions Indicates how many color- 
ant sets exist to render the given color. 

This number can be zero, one or it can be more than one. If the number is more than one, the set of colorants is selected 
25 that optimizes an additional criterion. According to a first embodiment, this criterion Is the amount of color change that 
results In the presence of colorant variations. According to a second embodiment, the one cdlorant solution Is retained 
that is connective to the colorant solutions for the colors along a selected path in color space. 

To invert the Neugebauer equations for an n-ink process with n>3. preferentially n-3 colorants are considered as param- 
eter and the n-ink model will be reduced to a Neugebauer process for three colorants that can be inverted by making 
20 use of the previous method. The choice of the colorant set to render a given color is obtained by minimizing a cost func- 
tion. 

The previous results can be extended for the localized Neugebauer model and hence by approximating a printer model 
with the localized Neugebauer model, it is possifc>le to Invert any n-ink process. 



25 Detailed description of the invention. 



The invention is described hereinafter by way of examples with reference to the accompanying figures wherein : 



Fig. 1 For a given X, Y and Z. each of the three Neugebauer equations represents a "surface" in the 3-dimensional 
30 Ci , C2. C3 space. The solution of the Neugebauer equations for a given X, Y and Z value can be seen as the 

intersection between the three surfaces that correspond to these values. In this example a cmy process is 
presented. The axis pointing to the right corresponds to yellow (y), the axis pointing to the left corresponds 
to magenta (m) and the third axis pointing upwards Is cyan (c). 
Fig. 2 The division of the colorant cube of a cmy process in 8 Neugebauer cells. 
35 Fig. 3 Newton Raphson iteration on one dimensional function. The iteration diverges because a point is hit in which 
the derivative has a zero value. 
Fig. 4 Example of a case where the three surfaces corresponding to the Neugebauer equations have more than 
one intersecting point. In this example surfaces in colorant space corresponding to constant tristimulus val- 
ues for a yellow, green and cyan ink process modeled by the Neuget^auer equations are represented. The 
«> axis pointing to the right corresponds to g, the axis pointing to the left corresponds to y and the third axis 

pointing upwards is c. 

Fig. 5 The section procedure of the colorant combination resulting in the most stable color rendering for a 3-lnk 
process modeled with the Neugebauer equations. 



45 While the present invention will hereinafter be described in connection with preferred embodiments thereof, it will be 
understood that it is not intended to limit the invention to those embodiments. On the contrary, it is Intended to cover all 
alternatives, modifications, and equivalents as may be Included within the spirit and scope of the invention as defined 
by the appending claims. 

In the following paragraphs, printer models are inverted for colors represented in a 3-dimensional space. 
so For up to three colorants, the inversion of the Neugebauer equations is presented. This solution can be extended for 
the localized Neugebauer equations and hence by approximating a printer model with the Neugebauer equations the 
Inversion of any printer model is obtained. In the case of multiple solutions, a selection mechanism Is presented. 
For four colorants, one colorant is considered as a parameter. 

This colorant is sampled between a minimum and maximum value and per sample point a 3-ink process is extracted 
55 from the 4-ink process by setting the colorant to its sampled value. To choose one colorant set. a given cost function is 
minimized. The methods to invert the Neugebauer equations, the localized Neugebauer equations and a general printer 
model are presented and several examples of cost functions are given. 

In the case of n colorants, n-3 colorants are considered as parameters, and per sample point in an n-3 dimensional col- 
orant space a 3-lnk process is inverted. Of all these colorant solutions, that solution is retained for which a given cost 



7 



NSDOCID: <EP 0763929Al„l_> 



EP 0 763 929 A1 



function is minimized. 
Inversion of a 3-ink process 
5 Neugebauer equations 

This solution is based on the observation that the set of Neugebauer equations for three colorants can be transformed 
into a 6th degree polynomial in one variable, and that robust numerical algorithms exist that, under atl conditions, find 
ail the roots (real and complex) of such a polynomial. The Neugebauer equations are here represented in the XYZ color 
10 space, but any other set of color coordinates in a 3-dimensional space can be used. The derivation of the 6*^ degree 
polynomial is presented here: By considering C3 as a parameter, the original Neugebauer equations can be rewritten 
as: 

X = (k o+k 3C 3) + (k 1+k 1 3C 3)0 1 + (k 2+k 23C 3)0 2 + (k 12+k 123C 3)0 , c 2 

15 

Y = (1 0+1 3C 3) + (1 1 +1 1 3C 3)0 2 + (1 2+1 23^3)^ 2 + (i 12+' 123^ 3)C 1 C 2 

M = (mo+m3C3) + (mi+mi3C3)C3 + (m2+m 3303)0 2 + ("112+"^ 123^3)^1^2 
20 Using the following definitions: 



Ko = 


ko + ksCa 


L-o = ^0 + '3C3 


Mo = 


mo + msCs 


Ki = 


ki + ki3C3 


Li = ti + I13C3 


h/li:= 


m-i + mi3C3 


K2 = 


^2 + ^*23^ 


L2 = I2 + 123^3 


M2 = 


m2 + m23C3 


Ki2 = 


= k^2 + *^123C3 


L^2 = '12 + ^23^3 


Mi2 = 


s m^2 + "^123^3 



30 the Neugebauer equations can be simplified to: 

X = Ko + K^Ci + K2C2 + K.,2CiC2 
Y = Lq + LiC^ + L2C2 + L,2C,C2 

35 

ZsMq + M^C^ -I-M2C2 + M^2^1^2 

This is a set of three equations with two variables. Assuming that ^ different from zero, the product c^C2 in the 
above equation can be eliminated: 

40 

(K o-X)M ,2 - (M o-Z)K 12 + c , (K ^ M ^g-M , K ,2) + C 2(K 2M la'M 2K ^2) = 0 
(Lo-Y)M,2-(Mo-Z)Li2-»-Ci(LiMi2-M,Li2) + C2(L2Mi2-M2Li2) = 0 
45 and, using the following definitions: 



Ao 


= {Ko-X)M,2 - (Mo-Z)Ki2 


Bo 


= (Lo-Y)Mi2-(Mo-Z)Li2 


Ai 






= (LiMi2-MiLi2) 


A2 


= (K2M12 - M2K12) 


B2 


= (^Mi2-M2Li2) 



the previous equation can be corrpactly rewritten as: 

55 

Aq -f c^A^ +C2A2 =0 
Bq + c,B^ +C2B2 «0 
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The solutions of these equations are immediately found using Cramer's rule: 



C2 = 



10 



By filling in and C2 in the third expression of the simplified Neugebauer equations, i.e. 

Z = Mq +M.,C^ + MgCg + M 12C1C2 
15 and working out the terms, the following 6**^ degree polynomial in C3 is okTtained : 



(M o-Z)det(D)(L , K 2-K 1 L 2)+(K o-X)det(D)(M 1 Lg-L 1 M 2)+(L o-Y)det(D)(K ^ M 

(Lo-Y)[(L2M ,2M 2)(K ,2-K i)+(K 12M2-K2M ,2)(Li2M i-LiM 

(L,K,2-L,2Ki)+(K2Li2-Ki2L2)(L,2M,-L,M,2)]+(Lo-Y)(Mo-2)[(K,2M2-K 

+(K2L ,2-K i2L2)(K nM ,2^ ,)]^ 0 



20 



25 



In which det(D) refers to the determinant of the matrix D : 



30 



^1 ^2 ^12 



35 



40 



This polynomial in C3 has up to six roots, some of which can be complex. They can all be found by means of. for exam- 
ple, the "Laguen-e" algorithm presented in William H. Press. Saul A. Teukdsky. William T. Vetterling, Brian P. Flannery, 
"Numerical Recipes in C". Cambridge University Press, 1992, Chapter 9.5. The corresponding values for Ci and <^ are 
found by filling in the roots C3 of the polynomial In the expressions 



_ ^2^0"^Q^2 



45 



SO 



55 



This leads to up to six triplets of c^ . Cg and C3 that meet the Neugebauer equations. 

If K|. Lj and are interchanged, exactly the same equations are obtained. From this follows that the roots of the 6th 
degree polynomial correspond to the solutions of the Neugebauer equations provided that K12 0 or L12 0 or M12 * 0. 
A special case occurs if a root would be found lor which K ,2 = L ^2 = 12 = ^ ■ ^^'^ would not necessarily be a 
solution of the original Neugebauer equations. In order to find out if this special case occurs, It Is evaluated if for one of 
the roots of the 6th degree K^2 = L,2 = M^2=0 ^"'y ^ is NOT the case, the root is also a solution of the Neuge- 
bauer equations. Othenwise the solution of the Neugebauer equations has to be searched by filling in the root value In 
the original Neugebauer equations, which In that case can be simplified to: 



X s k'o + k'^c^ + k'2C2 + k'i2C-|C2 
Y = I'o + 1' iC 1 + TgCg + r 12C1C2 
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Z = m^Q + m' -,c^ + m'gCg + 0^*12^1^2 

If one of the coefficients of the cross terms c^Cg is different form zero, suppose m'i2. the cross terms in c^Cg can be elim- 
inated in the first two equations. In this way a set of two linear equations are obtained in two variables, that can be 
5 solved using standard mathematical procedures. Only If for this CiC2-solution also the last equation is met. it is a solu- 
tion for the Neugebauer equations. A similar procedure can be applied if m'12 would be zero and one of the coefficients 
k*i2 or l'^2 »s different from zero. If all the coefficients of the cross terms C1C2 are zero, a set of linear equations are 
obtained that can be solved by making use of standard methods. 

It can be verified that all cases where at least one of the coefficients of the cross products between two colorants in the 
10 original Neugebauer equations are non zero, can be solved similarly by taking the remaining colorant as the parameter. 
In all other cases, the Neugebauer equations degenerate to a set of linear equations that can be solved using straight- 
foPA/ard techniques. The above explanation can be summarized as follows: 

the Neugebauer equations for a 3-ink process can be converted Into a 6*^ degree polynomial in one variable. Robust 
numerical methods exist to find the six roots of such a polynomial which can be real or complex conjugated. The roots 
15 of this polynomial correspond to six sets of three colorants with which a given color can be rendered if the Neugebauer 
equations are used as printer model. 

Hence the questions are answered how many solutions exist for the Neugebauer equations and what these solutions 
are. 

A Neugebauer process with less than three inks can be considered as a Neugebauer process with three colorants but 
20 a number of coefficients in the Neugebauer equations will be zero. Hence a Neugebauer process with less than three 
inks can also be inverted by making use of the inversion of a 3-ink Neugebauer process. 

Localized Neugebauer equations 

25 To invert the localized Neugebauer equations with up to three colorants, for each Neugebauer cell the equations are 
inverted and the real solutions that fall within the corresponding cell are retained. All these solutions are colorant sets 
that result in the required color. In most situations, a given color can only be obtained with one colorant combination. 
However, for some processes several solutions are found that belong to the same or different Neugebauer cells. 

30 General printer model 

Any printer model can be approximated by making use of the localized Neugebauer equations with any accuracy if a 
sufficient number of Neugebauer cells are used. Therefore it is possible to invert every printer model by approximating 
this model with the localized Neugebauer equations and inverting these equations. 

35 

Example 

By means of an example, it is demonstrated that the Newtonian iteration can sometimes "miss" the solution that is 
searched for. 

40 The following set of equations was obtained from measurements of a cmyk offset printing press. It predicts the XYZ 
color as a function of black (ci). cyan (03) and magenta (C3) colorants in combination with a fixed value of 100 % for the 
yellow colorant 

The coefficients are 'localized" and are only valid within the following ranges of the colorant values (in this paragraph 
the colorant values range from 0 to 1 .0): 

45 

0.2 < Ct < 0.4 
0.0 <C2 <0.1 
0.7 < C3 < 1 .0 

so 



55 
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59.72 


b = 


55.47 


mo = 


5.64 


ki = 


-26.74 


li = 


-37.10 


mi = 


-2.93 




-59.38 


»2 = 


-30.10 


tV\2 - 


14.22 


k3 = 


9.75 


'3 = 


5.06 


ITI3 = 


-16.69 


ki2 = 


-75.58 


Il2 = 


-66.93 


= 


-4.76 


ki3 = 


33.38 


Il3 = 


44.82 


mi3 = 


3.19 


k23 = 


73.89 


I23 = 


25.04 


m23 = 


-26.44 


kl23 


= 3.75 


h23 - 


14.06 


mi23 = 


30.94 



IS 

The equattonc pr edict that if the colorant values are set to 
Ct • 0 20 

Crj - 0 00 
20 Cj • t 00 

the resulting XYZ ccor is equal to: 

X •24S4 

25 Y.13% 
Z«240 

The mverscn at ihe Neugebauer equations by means of Newtonian iteration with starting values : c.|=c 2=0 3=0.5 
leads to the to^kwng soluton for the colorants : 

Ci « 0 38 

C2- -0 36 
C3 « 0 99 

35 By filling m these %«iijes m the equations it can be verified that although not physically realizable, this is indeed a valid 
solution from tne marhematical point of view. Inversion of the equations by means of polynomial root finding leads to the 
follawing four solutions (only the real solutions are retained): 



40 


Ci = 0.38 


Ci = 0.20 


Ci = 0.58 


Ci = 0.83 




C2 = -0.36 


C2 = 0.00 


C2 = 1 .06 


C2 = 0.04 




C3 = 0.99 


C3= 1.00 


C3 = 5.14 


C3 = 18.69 



45 

The first solution is Ihe one to wtiich the Newtonian iteration had converged to. Only the second solution falls within the 
colorant gamut in additon two other (non physically realizable) solutions were found. A pair of two complex conjugated 
solutions ts non chown 

50 Selection mechan»sm lor one set of colorant values 

To separate a color nto colorants, only one set of colorants is needed. Given the fact that the above method produces 
up to six solutions of ttw Neugebauer equations, a selection criterion is heeded to determine which solution is used for 
the separation of color mto colorant. 
55 In a first step the pav s of complex conjugated solutions that were found are rejected, since they do not correspond to 
physically realt^abie cokxant sets. In the case that all solutions are complex conjugated, the color lies in that part of the 
color space tfiat can *not t>e addressed" by the Neugebauer equations from colorant space. It is hence unprintable and 
for its rendering, a needs to be mapped onto a printable color. Examples of techniques to do so are described in US 
4,500,919, The Qarrut boundaries that are described in this patent can for example be calculated by means of the tech- 
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4 



niques described in the article by M. Inui "A Fast Algorithm for Computing Color Gamuts" (Color Research and Appli- 
cation. Vol. 18. No. 5. Oct. 1993). 

In a second step, these solutions of the Neugebauer equations are retained that lall inside the colorant gamut. For offset 
printing processes, e.g.. this range is usually from 0 to 100%. Assuming that the color is indeed a printable color, at least 
5 one such set of colorants exists with this property. If no such solution is found, the color Is unprintable and for rendering 
needs to be mapped onto a printable color. 

Till now. only solutions are retained that are both real and fall within the colorant gamut. In principle, any of these sets 
of colorants can be used for the correct rendering of the color, and the selection of one of them always contains an arbi- 
trary element. 

10 Two approaches are presented here but many more exist and iaW within the spirit and the scope of the Invention. 

According to a first approach that solution is selected that yields the most stable rendering of the color on the printer. 
With most stable rendering is meant that set of colorants that in the presence of incremental changes on the colorant 
values yields the smallest change In color. The change of color as a result of colorant changes is easily calculated by 
evaluating the total derivative of the forward Neugebauer model for the particular colorant set. i.e. 

75 

dX/dc^ dx/dc^ dX/dc^Y^^ ^ 
dY/dc^ dY/dc^ dY/dc^ dc^ 



20 



dx 

dY 

dz 



25 By calculating the singular value deconrposition (SVD) of the previous 3X3 matrix, an orthogonal axis transformation is 
performed in such a way that one of the new axes lies along the maximum color change for a unit change in the colorant 
space. Tliat color change corresponds to the largest singular value, because it is expressed in color units per unit col- 
orant change (the normal Euclidean distance measure is used in both colorant and color space). The colorant combi- 
nation for which the largest singular value is the smallest is retained. The selection of the colorant combination resulting 

30 in the most stable color rendering is presented in the flowchart in Figure 5. The SVD of a matrix can be calculated using 
the method presented in William H, Press, Saul A. Teukolsky, William T Vetlerling. Brian R Flannery, "Numerical Reci- 
pes in C". Cambridge University press, 1992. 

According to a second approach, continuity of the colorant changes is given preference along certain paths in color 
space. Especially in areas of smooth tone and color transitions, abrupt spatial changes in colorant values are highly 

35 undesirable since they create visible artifacts in the presence of misregistration between the colorant planes. For a cmy 
printing process, for example, this can be translated into the requirement to make the rendering of degrades from white 
to cm. cy and my and from cmy to cm. cy and my as smooth as possible. In the subsequent calculation of the colorant 
values for the rendering of such color paths, always the set of colorants that differs the least (In terms of least mean 
square) from the set of colorants at the previous point on the path Is selected. The same principle can be used when 

40 the points of a 3-dlmensional separation look up table are to be filled in. 

Inversion of a 4-ink process 
Neugebauer equations 

45 

Because the Neugebauer equations are continuous functions, there may be up to six paths in colorant space that result 
in the same color. This can be seen as follows. If the fourth colorant Is kept constant, there may be up to six solutions 
to render a given color. By varying the fourth color inf initesimally. the previous six solutions will also be changed inf ini- 
tesimally due to the continuous relations between the color values and colorants. Hence by changing the fourth color- 
so ant. six paths are obtained in the colorant space that result In a given color. 

To select one colorant combination, the minimum and maximum value for one colorant is determined with which a given 
color can be achieved. These values can for example be calculated by inverting all the extracted lx)undary 3-ink proc- 
esses of the 4-lnk process for the given color. Per colorant, the minimum and maximum value of all these colorant solu- 
tions Is the minimum and maximum value of this colorant with which the given color can be reached. The six paths can 
55 now be sampled by choosing a value for one colorant between its minimum and maximum value, extracting the 3-ink 
process by filling put this value in the 4-ink process, and inverting the extracted 3-ink process for the given color. From 
all these colorant combinations that fall inside the colorant gamut the one set is retained for which a given cost function 
is minimum. Several examples of cost functions are: 
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Most stable rendering : In the case of a 4-ink process the relation between color changes and colorant changes is 
expressed with a 3X4 matrix. The colorant set for which the largest singular value is minimum is retained. 
Minimum colorant : the colorant combination for which a given colorant value is minimum 
Maximum colorant : the colorant combination for which a given colorant value Is maximum 
5 • Percentage colorarrt the colorant combination tor which a given color corresponds with a percentage between its 
minimum and maximum value 

Minimum colorant sum : the colorant combination for which the sum of the colorants is minimum 
Maximum colorant sum ; the colorant combination for which the sum of the colorants is maximum 
Equal colorants : the colorant combination for which the difference between the colorant values is as small as pos- 
10 sible 

Moire : the colorant combination that is the least sensitive to moire 

Path : choosing one path out of the six that results in the given color and applying a cost function to select one col- 
orant combination on this path. 

Continuity : colorant combinations are chosen in such a way that continuous color changes are separated into con- 
15 tinuous colorant changes. 

It is dear that this list of cost functions is unlimited and that many more cost functions exist that fall within the scope and 
the spirK of this invention. 

20 Localized Neugebauer equations 

The previous method can be extended for the localized Neugebauer equations with four colorants. Per Neugebauer cell 
there may be up to 6 paths that map to the same color in color space and hence there may be several paths in the col- 
orant domain that result in the given color if all paths in the Neugebauer cells are taken together. However, because the 
25 localized Neugebauer equations are continuous functions, there will be paths of neighboring cells that are connected 
with each other. 

To calculate the minimum and maximum value of a colorant to obtain a given color, a distinction has to be made about 
the behavior of the printer. For well-behaved processes, such as a cmyk offset printing press, the minimum and maxi- 
mum value of a colorant to form a given color can be found by determining the extracted boundary 3-ink processes of 

30 the localized Neugebauer and inverting these extracted boundary 3-ink processes for the given color. Per colorant, the 
minimum and maximum value of these colorant solutions are the minimum and maximum value tor this colorant with 
which the given color can be reached. For some unconventional processes however, previous procedure has to be 
applied per Neugebauer cell. i.e. determining the extracted boundary 3-ink processes and inverting these extracted 
boundary processes for the given color Per colorant, the minimum and maximum colorant value of all the previous col- 

35 orant solutions are the minimum and maximum value of these colorants with which the given color can be obtained. 
Then, one colorant is sampled between its minimum and maximum value, the 3-ink process is extracted by filling out 
this value in the 4-ink process, and finally the extracted 3-ink process is inverted for the given color. From all these col- 
orant combinations that set is retained that falls inside the colorant gamut of the n-ink process and tor which a given 
cost function is nninimum. 

40 

General printer model 

Any printer model can be approximated by making use of the localized Neugebauer equations with any accuracy if a 
sufficient number of Neugebauer cells is used. Therefore it is possible to invert every printer model by approxtnnating 
45 this model with the localized Neugebauer equations and inverting these equations. 

Inversion of an n-ink process 

Neugebauer equations 

so 

Previous methods can be generalized for any number of inks. The minimum and maximum colorant value to reproduce 
a given color is obtained by determining ail the extracted boundary 3-ink processes of the n-ink process and inverting 
these extracted boundary 3-ink processes for the given color. Per colorant, the minimum and maximum value of all 
these colorant solutions are the minimum and maximum colorant value for this colorant to reproduce the given cotor. 
55 For an n-ink process with n > 3 there are n-3 degrees of freedom, and hence n-3 colorants will be considered as param- 
eters, that can have colorant values between their minimum and maximum value with which a given color can be 
reached. This (n-3)-dimensional parameter space is sampled and the 3-ink process is extracted form the n-ink process 
by setting the (n-3) colorants in the printer model of the n-ink process to their sampling values. This extracted 3-ink proc- 
ess is inverted and for all the sampling points that solution is retained for which a given cost function is minimum, taking 
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into account the colorant iimitations of the colorant gamut of the n-ink process. 
Localized Neugebauer equations 

The previous method can be extended for the locaftzed Neugebauer equations with n colorants. 
To caicUate the minimum and maximum value of a colorant to obtain a given color, a distinction has to be made about 
the behavior of the printer. For well-behaved processes, such as a cmyk offset printing press, the minimum and maxi- 
mum vai je of a colorant to form a given color can be found by determining the extracted boundary 3-ink processes of 
the iocalzed Neugebauer and inverting these extracted boundary 3-ink processes for the given color. Per colorant, the 
minimum and maxirrujm value of these colorant solutions are the minimum and maximum value for this colorant with 
w^lc^ the given color can be reached. For some unconventional processes however, previous procedure has to be 
applied per Neugebauer cell. i.e. determining the extracted boundary 3-ink processes and inverting these processes for 
the given coky Per colorant, the minimum and maximum colorant value of all the previous cok^rant solutions are the 
mmirrxim and maximum value of these colorants with which the given color can be obtained. 

Then, n-3 cok>ants are sampled between their minimum and maximum value, the 3-ink process is extracted by filling 
out tfies^t? values m\ the pr inter model of the n-ink process, and finally the extracted 3-ink process is inverted for the given 
color From aU ihese cokxant combinations, that set is retained that falls inside the colorant gamut of the n-ink process 
arxJ f r which a gt^ers cost function is minimum. 

General printe* rnodo* 

Any printer mocM can t% approximated by making use of the localized Neugebauer equations with any accuracy if a 
suff k:ient numtMr o# Neugebauer cells is used. Therefore it Is possible to invert every printer model by approximating 
this model with th# tocakzed Neugebauer equations and inverting these equations. 

Having described wy <U4Aii preferred embodiments of the cun-ent invention, it will now be apparent to those skilled in the 
art that numerous modrfications can be made therein without departing from the scope of the Invention as defined in 
the following claini& 

Claims 

1 • Method of oomputng a set of colorant values for obtaining a given color by a color reproduction device comprising 
the iolk>wing steps 

modeling said cokx reproduction device by Neugebauer equations ; 

f inckng a piurab?y of sets of colorant values by inversion of the Neugebauer equations for said color ; 
seiectng one set of cok>rant values from sakj plurality of sets. 

2. Method accorc^ to daim 1 . further comprising the following steps : 

def ning a colorant space ; 

dtvding ssud colorant space in Neugebauer cells : 

modeliing sax) cokx reproduction device by a set of Neugebauer equations per Neugebauer cell ; 
finding a gkMakry of sets of colorant values by inversion of the Neugebauer equations for said color ; 

- seiectng or^ set d colorant values from said plurality of sets, said colorant values being comprised within the 
corresporxlir^ Neugebauer cell. 

3. Method accocdng to any of daims 1 or 2, further comprising the following steps : 

extracting a 3-4nk process from said Neuget>auer equations : 
transfornwig sa«1 3-ink process to at least one polynomial in one variable ; 
inverting sa«l 3-ink process by finding a plurality of roots of said polynomial. 

4. Method acoordng to claim 3. wherein said 3-ink process is extracted from an N-ink process by taking N-3 colorants 
as parameter 

5. Method according to daim 4. further comprising the following steps : 

- finding ky at least one parameter a minimum and maximum colorant value, suitable for reproducing saki given 
color : 
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- sampling said parameter within a range between said minimum and said maximum. 

6. Method according to any of claims 1 to 5. wherein said set Is selected by minimisation of a cost function. 
5 7. Method according to any of claims 1 to 6. further comprising the following steps : 

rejecting complex or imaginary solutions ; 

- rejecting solutions outside an established colorant space. 

10 8. Method according to any of claims 1 to 7, wherein unprintable colors are mapped to printable colors within a color 
gamut of said color reproduction device. 

9. Method according to any of claims 1 to 8. further comprising the following steps : 

fs - assessing for each set off colorant values a color change caused by Incremental changes of colorant values : 
selecting the set of colorant values with the smallest color change. 

10. Method according to any of claims 1 to 8. further comprising the following steps : 

cc selecting at least one continuous path in color space ; 

- selecting the set of colorant values situated on a continuous path in colorant space, that is transformed to said 
path in color space by the Neugebauer equations. 

11. Method according to any of claims 1 to 10. wherein said Neugebauer equations are trilinear polynomials. 

12. Metfxxl according to any of claims 2 to 11 . wherein said Neugebauer cells are mutually disjunctive and preferen- 
tially cuboid. 
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